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Abstract 



We consider 3D active plane rotators, where the interaction be- 
C^ I tween the spins is of XY-type and where each spin is driven to rotate. 

For the clock-model, when the spins take A^ ^ 1 possible values, 
we conjecture that there are two low-temperature regimes. At very 
low temperatures and for small enough drift the phase diagram is a 
small perturbation of the equilibrium case. At larger temperatures the 
massless modes appear and the spins start to rotate synchronously for 
arbitrary small drift. For the driven XY- model we prove that there is 
essentially a unique translation-invariant and stationary distribution 
despite the fact that the dynamics is not ergodic. 
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1 Introduction 

Understanding nonequilibrium phase transitions is a major challenge of sta- 
tistical physics, bearing many different aspects. Today not much of a sys- 
tematic theory exists, with few experimental and even fewer mathematical 
results. One important question of this multifaceted subject is to describe 
the changes to the equilibrium phase diagram when a steady nonequilibrium 
driving is added. The best known examples are driven diffusive lattice gases 
- like boundary driven and asymmetric exclusion processes, see e.g. the 
books [m H] , including many results of computer simulations. The nonequi- 
librium there originates from installing differences in chemical potentials or 
from adding nonconservative external fields, and follows the prescription of 
local detailed balance, [15]. In the present paper a uniform nonequilibrium 
force drives the internal degrees of freedom, the planar rotating spins. The 
spins take values on the unit circle and are placed on the sites of a regular 
lattice. When mutually uncoupled, all spins undergo the same non-reversible 
Markov evolution with a bias in the direction of rotation. The interaction 
couples nearest neighbors, x ~ y, with energy following the XY-model, 

H{<^) = -^cos{<^^-<^y) (1) 

for "angles" (px,x & iP . The possible values of the angles determine the 
nature of the spins. A first choice is to take '^^ = 2'Kk/N, k = 1,2, . . . ,N on 
the discrete circle with A^ possible values. For N = 2 that is the Ising model; 
for A^ = 3 it is equivalent with the g = 3 Potts model. The second possible 
choice is formally obtained in the A^ ^ +oo limit, and has a continuum of 
values (px G [0, 27r]. That truly corresponds to the XY-model where the spins 
are plane rotators having unit length. 

The purpose of the paper is to discuss the modification of the phase diagram 
when a nonequilibrium driving is inserted that induces biased rotation of the 
spins over the circle. By doing so the following phenomena can be addressed: 

a) the uniqueness of the stationary distribution accompanied by breakdown 
of ergodicity - in the sense that some initial data do not relax; 

b) the presence of macroscopic dynamical coherence; 

c) the stability of equilibrium phases against small nonequilibrium driving; 
We first briefly introduce each of these points, to realize them more concretely 
in later sections. 



1.1 Unique stationary distribution without ergodicity 

Are there stochastic dynamics with a unique stationary distribution, which 
are not ergodic? Without any further restrictions this question is easy and 
not very interesting, with the answer being 'yes'. For discrete time a simple 
example is given by the two-state Markov chain 

+ 1^-1^+1^..., (2) 

flipping deterministically at each time. Similarly, in continuous time we can 
take the rotation over the circle with a constant angular speed v. 

e{t) = e{Q)+vt mod27r (3) 

One wonders whether one can construct non-degenerate random processes, 
which exhibit the above prototypical behavior. Of course, we must then 
consider infinite- volume interacting particle systems and infinite probabilistic 
cellular automata, since finite state non-degenerate Markov processes are 
always ergodic. 

For the case of probabilistic cellular automata (discrete time parallel up- 
dating of spins), a construction, mimicking the behavior ([2]) is presented in 
a recent paper by [3] . However, the example of [3] still has some degeneracy, 
because for every time T one can present two local events, A and B, such 
that the transition probability px (^l-B) in T steps vanishes. Thus, we feel 
that a truly non-degenerate discrete time example is still missing. We be- 
lieve that the discrete time version of our 3D driven clock model gives such 
non-degenerate example. 

Our constructions below present the case of rotating interacting spins in 
continuous time. The rotation speed f in ([3]) will be induced by the nonequi- 
librium driving, and the angle 9 should be thought of as the order parameter, 
or collective phase, of the model. The dynamics will be nondegenerate, as 
local fiuctuations in the phase are allowed. See Theorem [7] for the precise 
result. 



1.2 Macroscopic coherence 

Not surprisingly, the mechanism above connects with the old but still not 
completely resolved question of whether macroscopic dynamical coherence 
or pattern formation in spatially extended systems can be obtained by local 



translation-invariant and non-degenerate updating of spins and whether that 
is even possible for a continuous time (sequential) dynamics and for dynamics 
that satisfy detailed balance, [H |5]. In fact, an example with that flavor was 
recently described in |18] . There an infinite queuing network was considered, 
with several types of clients and with exponential service times. In the high 
load regime the system exhibits coherent behavior. That means that if the 
initial state of the network is close to the 'coherent' one, characterized by 
a given value of the 'phase' observable (which takes values on the circle), 
then in the process of evolution this phase evolves with a constant speed, is 
never 'forgotten', and the initial synchronization is never broken. Still the 
system has a unique stationary distribution, with the phase being uniformly 
distributed over the circle. In the language of queuing networks it is an 
example of violation of the Poisson hypothesis. Yet, this example lives not 
on a lattice with short range interactions, but on a mean-field graph (which 
is, in some sense, an infinite complete graph). 

There is also a vast literature on the emergence of synchronized rotators using 
variants of the so called Kuramoto model; for a review, see pQ. Recently, a 
mean field analysis for active rotator models was carried out in [12] , which for 
some choice of the drift is the mean field version of the model we consider later 
in ([7]). We believe that rotating states emerge in low temperature uniformly 
driven A^-clock models, if A^ is sufficiently large. See the conjectures in 
Section [3J 

1.3 Stability of equilibrium phases 

When the equilibrium model has finitely many macroscopic phases in some 
regime of its parameters, then we expect that these remain in place for small 
driving. Below a critical driving, the changes in basin of attraction and in 
macroscopic appearance will be small. To understand the nature of that crit- 
ical driving, one must realize that the stability of equilibrium phases requires 
some sort of free energy barriers, or, in particle language, the excitations must 
be massive. Therefore, the presence of soft modes, or Goldstone bosons, can 
break stability. In the context above, that means that the critical driving 
(the minimum we need to truly disturb the phase diagram) will go down as 
N '\ +CXD. See conjectures |3]-|1] for more precise speculations. 

The following Section [2] contains the details of the nonequilibrium model 
- driven A^-clock models - together with a summary of the situation in equi- 



libriuin. Section E] is devoted to what we believe happens for finite A^; these 
are mostly a collection of conjectures in which we firmly believe but where 
the proofs are missing. Some of this is remedied in Section H] for the driven 
XY-model where the picture is more complete. The main result is that the 
3D driven XY-model shows nonergodicity, while having a unique stationary 
translation invariant distribution, at (almost) all low temperatures. 

2 The iV-Clock Model 

The A^-clock model is an interacting particle system that lives on T? . At each 
site X G Z^ there is a spin a^^ G Zat, where Z^r C §^ C C^ is the group of A^-th 
roots of unity. Each spin Gx has its clock, and when the clock rings, the spin 
jumps to one of the two 'nearest' values: Ua; — )■ a^ = exp {±^} a^. That 
is equivalent to introducing the angles i^^. with o^ = expiip^, ip^ = 2'7ik/N, 
k = 1,2, ... ,N, and moves y^^. — )■ yja- ± 27t/N. In what follows we will use 
the notation ( (</?) for exp {iip} , ( (ip) G S-*^ C C^, so in particular 

a^ -^ a^ = C i± — j a^. 

The particles are interacting, with the energy given by ([T]). We define the 
rates c (x, a, ±) of the jumps a^ — ?• cf^ of the spin a^ in the environment a 
at inverse temperature /3 by 



c{x,a,±) =p±exp«{ -^ 



CO?.[iPx- Lpy±—] - COS (V^x - ify) 



(4) 



(sum over nearest neighbors y of x G Z^), where the numbers p+ > P- > are 
two extra parameters. Their difference is measured hy d := log (p+/p_) > 
and is called the drift. One can imagine it as the coupling of the planar 
rotator with a magnetic field that acts perpendicular to the plane. We call 
the above model the W-Clock model' with a drift. (Of course, the drift does 
not make sense for N = 2 (Ising model).) 

We first describe the properties of the symmetric Clock model, when 
the drift d = 0. Note that in this case the evolution defined above satisfies 
detailed balance. That is what we call the equilibrium or symmetric Clock 
model. Then, the Gibbs measures for ([1]) are reversible stationary measures. 
We can thus use the results of [5], theorem 4.6: 



Theorem 1 (symmetric Clock model) There exists a value /3o of the inverse 
temperature, such that for every (3 > [3q and every N > 2, the symmetric N- 
Clock model has at least N different extremal stationary distributions, (■)^^ ^ , 
(^fc = C (^) 5 ^ = 1,...,A^. These states are translation-invariant, exhibit 
long-range order and are magnetized: 

^^^hk,P = "^^ (^) ^^^' ^^^^ ^^ *^^) > for /3 > ISq. 
(Here we interpret the spins ax as elements of C^ .) 

Observe that the finite Pq above remains the same for all A^ > 2, which 
makes the Theorem unreachable for the standard low-temperature analysis 
based on the Peierls condition. For example, the Pirogov-Sinai theory pn[2T] 
would establish the stability of the A^ ground states only for /3 > (3^^ , where 
/3^^ -> oo as A^ — )■ oo. The reason for that is not purely technical: indeed, 
in the domain of validity of the Pirogov-Sinai theory one necessarily has 
additional properties of the pure phases, such as the exponential decay of 
the truncated correlation functions. However, we believe that in reality such 
exponential decay holds only for low enough temperatures, (3 > (3^, with 
/3^ — )■ oo as A^ — )■ oo, so the PS-method can not be improved to reach (3q. 
Moreover, we think that the following is true: 

Conjecture 2 There exists a value N, such that for each N > N the 3D 
symmetric N -clock model undergoes two phase transitions. Namely, for all 
(3 > (3^ it has N pure magnetized phases, with exponential decay of truncated 

5^ 



correlations, with /3^ — )■ oo as N ^ oo. For smaller intermediate values of (3, 



(3^ > f3 > /3^ it also has at least N pure phases with non-zero magnetization; 
however, the correlation decay in these phases is only algebraic. (A stronger 
recent conjecture of JB^ even talks about a continuum of pure phases, {■)cri 
C G S-'^.y) Finally, for (3 < f3'^ the model has only one Gibbs state, again with 
exponential decay of correlations. 

We know from [2] that there is no intermediate phase for the Ising model, 
hence N > 2. The conjectured behavior is somewhat similar to the one 
for 2D Clock models; it was proven in [TT] that these indeed undergo a 
Berezinskii-Kosterlitz-Thouless phase transition. There, massless Goldstone 
modes appear (hence our notation /3^). 



3 Conjectures for the driven Clock-model 

We now discuss the situation with non-zero drift (i, (where we do not have 
detailed balance). Let us run our A^-Clock model with a drift for a time 
duration T, starting in one of the equilibrium phases (■)cfe,/3. Let us denote 
the resulting state by (•)o'^ . 

Conjecture 3 For every N > N, P > /3^ there exists a critical value 
dcr (/3, N) of the drift, < d^r [P, N) , such that the following holds: 

1. if \d\ < dcr{f3,N), then the state (■)^'^ approaches the state {■)pd, as 
T — )• oo, which is magnetized: 

(^o)^,, 7^ 

and is close to (■)(jj.,/3 for small d, in the sense of expectations of local 
observables; 



k,T 



m 



2. if d > dcr (/3, N) , then the state (■)a'(^ is a 'rotating' state as T —)■ oo (• 
particular, it has no limit as T ^^ oo/ Namely, there exist two periodic 
functions: m(T) = m (T; /3, iV, rf) > and $ (T) = $ (T; /3, N, d) , z. e. 

m{T + uj)=m{T), $ (T + w) = $ (T) , 

with period u being the mean angular velocity, u = u {/3, N, d) , and a 
phase shift 0^ = 0^ (/3, A^, d) , such that 

(^x) J'J - m (T) e^(*(^)+'^^) ^ as T ^ oo. (5) 

(Here we again are treating the spin ax as belonging to C^.) 

It is interesting to compare the curve of the states (■)«'^ , T > 0, with the 
conjectured ([6]) pure phases of the symmetric clock-model, (■)r/3! C ^ S^- 
One cannot be stopped from guessing that perhaps for every T we have 
{■)m = (■)c(T),/3 ' for some C (T) = C (T, k, d, /3) G S^. 

The next conjecture deals with the behavior of the critical drift dcr intro- 
duced above. 

Conjecture 4 The critical drift is positive at low temperatures. It decreases 
to become zero at (3^: 

G 



dcriP,N) = OforP'^<P<P'^. 



The rationale behind this conjecture is that at temperatures above (/3^) 
the A^-Clock model enters into the spin-wave phase or Goldstone modes 
regime and so qualitatively should behave like the XF-model, which is in 
the rotating phase for any non-zero value of drift, see below. This similar- 
ity of the intermediate phases with the XY-raodel is the basis of all our 
speculations. Hence, we believe in the following 

Conjecture 5 For every /3 large enough there exists N = N (/3) , such that 
for any N > N (/3) and for all d > the N- Clock model with a drift d has a 
continuum of different rotating states. 

One might wonder whether there is a difference between the structure of 
the stationary states in the Pirogov-Sinai regime and in the Goldstone modes 
regime, when d > d^r (/?, N) , i.e., when we are in the regime of rotating states. 
We expect the answer to be positive: 

Conjecture 6 1. Rotating soft modes. In the regime /3 G (/3^,/3^), 
d> there is a unique stationary distribution, (■)o^. It is translation 
invariant and has zero magnetization. It is given by the limit 

<-)« = .l-/r"<-)S''r (6) 

(which does not depend on k). 

2. Rotating PS. In the regime f3 > /3^, d > dcr{P,N), in addition to 
the time-stationary translation invariant state (0) there are also time- 
stationary non-translation invariant states (the 'Dobrushin' states). They 
are given by the same formula i^ , where instead of the states (■)a'^ 

one should use the states (■)^^' . The latter are obtained by starting 
the driven N- Clock dynamics with the measure 6±k that gives weight 1 
to the configuration 

( g2ufci/7V j^^ Xi > 

and where we suppose for simplicity that N is even. 



The Dobrushin time-stationary non-translation invariant states have a 
rigid interface at the level Xi = 0. In a typical configuration drawn from such 
a state the spins on different sides of the interface are pointing in (approxi- 
mately) opposite directions, though the direction itself can be arbitrary. One 
should remember here that no Dobrushin states exist in the 3D XY-model, 
as shown in [TD]; see also the discussion in [20] ■ In fact, it is argued in [TU] 
that there are no non-translation invariant states at all in the 3D XY-model. 
This is the basis of our Conjecture |6l More precisely, we believe that in 
all cases there is a unique translation invariant stationary distribution for 
/3>/3]v^c^>4r•(/3,iV). 

4 The three dimensional Xy-model: main re- 
sult 

The dynamical Xy-model with a drift - called dXY model below - can 
be obtained from the energy ([T]) , and the dynamics (jl]) by taking the limit 
A^ — )■ cx), in a diffusive rescaling of time by t — )■ t/N^. Alternatively, it is a 
3D model of coupled Brownian motions ipx,x G Z^, on circles ipx G S^ C C 
The Brownian motions (fx have a constant drift, d, and they are interacting 
via the nearest neighbor attraction ([I]) . The dynamics for ipx{t) G [0, 27r] is 
then as follows: modulo 27r, 

d^x{t) = ddt-j—dt+J- dWx{t), (7) 

dipx V P 

where the Wx{t) are independent standard Wiener processes and 

dH 



Y^ sm{ipx{t) - ipy{t)) 



y:y~x 

The formal generator of this process, acting on local smooth functions / is 

df 



Lf = Lof + dJ2 



d^x 



i„;.53[_|^|ia|!(] (s) 



Observe that Lq commutes with the new generator for the driven modeh 
[L, Lq] = because for all x, 



E 



S'H ^ 



difxdify 



Of course L — Lq generates independent rotation on each angle with an- 
gular speed d and thus commutes with the generator Lq of the (undriven) 
XF— model. In other words, the dynamics of the XF— model can be inter- 
changed with uniform rotation of all spins. 

Much of the equilibrium structure of the XF-model is known. At low 
temperatures /3~^ the 3D XF-model has a continuum of translation-invariant 
Gibbs states. They can be obtained as thermodynamic limits {■)ra of the 
finite-volume Gibbs states with coherent boundary conditions ip^ = (^ & 
§^. outside the volume. These translation invariant states have non-zero 
spontaneous magnetization, 

(<^o)c,/3 = "^ (/3) C, with m (/3) > 0, 

see [g. 

Here is our main result. Consider the set S of stationary and translation- 
invariant distributions for the dXY— model. 

Theorem 7 The set S is a singleton for almost all temperatures, while at 
sufficiently low temperatures there exist rotating states as in Conjecture {^2: 

(■)/3',d = (■)c+dT,/3- 



Proof. We start by repeating that the evolution of the random variables 
i^x (t) = fx{t)—dt is that of the XF- model with zero drift. If therefore /i is a 
stationary distribution for the dXY— model, then fj, is periodically repeated 
under the symmetric Xy-dynamics. Suppose now that /i is translation- 
invariant. Then, Holley's argument shows that /i is in fact a translation 
invariant Gibbs measure for the XY- model, [13]. Moreover, fi must then be 
rotation invariant (S^-invariant) since stationary states of the dXY —model 
have zero magnetization. From pj6] it then follows that fi is unique for almost 
all temperatures. 

On the other hand, at low temperatures /3~^ the 3D XF-model has a con- 
tinuum of translation-invariant Gibbs states. These phases of the XF-model 
correspond in an evident way to rotating states of the dXY—m.odel. These 

10 



rotating states then do not converge to a stationary state of the dXY modeL 

■ 

Remarks: 



1. The Gibbs field (■)^ , defined by 



/3 



(■)«= {■)cBdC 



is §^-invariant. We beheve that for all /3 the state (■) ^ is the only translation- 
invariant Gibbs state of the Xy-model, which is §^-invariant, which would 
remove the "almost all" . (This statement is proven to hold for the XY-raodel 
for almost all values of /3, |1UJ.) 

2. Note, however, that by changing the interaction from cos {ipx — ^y) to 
/ — cos{i^x-'^y) \ _ J g^ i^y passing to the so-called 'very-nonlinear a— model' - 

we obtain an example of a system which at some temperatures has at least 
two translation-invariant S^-invariant Gibbs states (for p large enough), see 

m- 
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